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The phenomenological symplectic model with a Davidson potential is used to construct rotational 

states for a rare-earth nucleus with microscopic wave functions. The energy levels and E2 transitions 

obtained are in remarkably close agreement (to within a few percent) with those of the rotor model 

with vibrational shape fluctations that are adiabaticalUy decoupled from the rotational degrees of 

freedom. An analysis of the states in terms of their SU(3) content shows that SU(3) is a very poor 

dynamical symmetry but an excellent quasi-dynamical symmetry for the model. It is argued that 

^^ I such quasi-dynamical symmetry can be expected for any Hamiltonian that reproduces the observed 

^\ . low-energy properties of a well-deformed nucleus, whenever the latter are well-described by the 

0^ ' nuclear rotor model. 
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I. INTRODUCTION 



A microscopic theory of nuclear structure would be very incomplete without a satisfactory description of nuclear 
^ I rotational states in terms of many-nucleon quantum mechanics. However, while the states of a truly rigid rotor can 
O^ ' be handled with ease [|lj , they do not have square- integrable wave functions in either a spherical vibrational- model or 
^«J ' many-nucleon Hilbert space. Moreover, the expansion of liquid-like, soft-rotor, wave functions on any spherical basis 
if-^ : is slowly convergent. This means that many major shells are required for a realistic shell-model theory of nuclear 
f^ ' rotational states. It also means that a realistic calculation of nuclear rotational states in terms of interacting nucleons, 
^s^ , without a priori knowledge of the kinds of correlations to expect, is an impossibly difficult task. The fact remains 
0^ that nuclear rotational bands are exceedingly simple; they are essentially characterized by a few intrinsic quadrupole 
(-H [ moments and moments of inertia. Furthermore, they are observed in a wide range of nuclei throughout the periodic 
table. Thus, it would appear that the phenomenon is remarkably robust and relies very little, for its existence, on the 
details of the two-nucleon interaction. 

As demonstrated by calculations within the framework of the nuclear symplectic model [^]-Q, it is possible to 
construct model rotational states with many-nucleon wave functions. The problem is to understand why the model 
works as well as it does; i.e, why the many residual interactions that strongly break symplectic symmetry, by mixing 
states of different irreps, do not destroy the predictions of the model. 

To understand this, the effects of the dominant symmetry breaking interactions have been explored one at a time. 
; I It has been shown that, while a spin-orbit interaction may mix SU(3) irreps strongly, it does so, in large-dimensional 

C^ ^ irreps, in a highly coherent way such that the rotational spectrum of the model survives g. Similar coherent mixings 
of SU(3) irreps by pairing forces were found while investigating the transition of a many-fermion model M from a 
superconducting phase to a rotational phase with variation of the relative strengths of short-range (pairing) and long- 
range (quadrupole-quadrupole) interactions M. In this paper we examine in some detail the very strong mixing of 
major harmonic-oscillator shells and SU(3) irreps, within the framework of the symplectic model. Again the mixings 
are extraordinarily coherent. We refer to this coherent mixing as a quasi- dynamical symmetry WM- 

We conjecture that, because the residual interactions separately preserve quasi-dynamical symmetry, they continue 
to do so when combined for suitable ranges of their strengths. This conjecture is strongly supported by the observation 
of rotational states in nuclei which demonstrates unequivocally that they do survive. 

The discovery of quasi-dynmical symmetry gives optimism that realistic calculations of rotational states, which 
simultaneously take into account the mixings between major shells and, for example, the effects of spin-orbit and 
short-range (e.g., pairing) interactions, may be possible. Up to the present time such calculation have not been 
possible except within the framework of a highly constrained (e.g., Hartree-Bogolyubov) approximation. 

The calculations reported here are carried out with a simple Davidson interaction potential [g|. This interaction 
preserves a higher, symplectic model, symmetry which makes calculations in a very large multi-shell space possible. 
With this interaction, the Hamiltonian for a diatomic molecule was recently diagonalized by a simple non-linear 
transformation from a spherical vibrator basis to a soft rotor basis using an su(l, 1) spectrum generating algebra 



(isomorphic to sp(l, M)) [|T^. A similar transformation ||l^,|l2l gives rotational states in the Bohr-Mottelson collective 
model [ij]. 

For a many-nucleon nucleus, one can construct a Hamiltonian comprising the many-nucleon kinetic energy, a spher- 
ical shell model potential, and the same Davidson potential as used to obtain rotational states in the pheomenological 
Bohr-Mottelson model. The Davidson potential is a scalar function of the nuclear quadrupole moments expressed in 
terms of nucleon coordinates. Thus, it is microscopic and rotationally- invariant. More importantly, it is expressible in 
terms of an sp(3,M) spectrum generating algebra and and gives a Hamiltonian that is diagonalizable within a single 
irrep of the symplectic model f^-^- We are not able to give analytical expressions for its wave functions, as we were 
pU for the diatomic molecule and Bohr-Mottelson model. But, we are able to compute its states numerically and 
expand them on a spherical shell-model basis. This is done by diagonalizing the Hamiltonian in a basis that reduces 
the dynamical subalgebra chain 

sp(3,M) Dsu(3) Dso(3), (1) 

where su(3) is the symmetry algebra of the spherical harmonic oscillator shell model |Q . 

Symplectic model calculations in an su(3) D so (3) basis have been developed to a fine art using VCS (vector 
coherent state) methods [ |l4| , |l5[ to calculate matrix elements. Furthermore, they can be carried out in large spaces 
with relative ease. 

An important property of the sp(3,R) algebra is that it contains, as a subalgebra, a spectrum generating algebra 
for a rigid-rotor model, namely Ui's [E^]so(3) algebra mn. However, one cannot diagonalize a symplectic model 
Hamiltonian in a basis that reduces the subalgebra chain 

sp(3,M)D [M^]so(3)dso(3), (2) 

because, as mentioned above, rigid-rotor wave functions are not expandable in a shell-model basis (in technical 
terms, the decomposition is not a direct sum; it is a direct integral.). One notes, however, that the Hamiltonian 
we are considering as a model of a soft rotor has two components: a spherical shell model part which reduces the 
sp(3,R) D su(3) D so(3) chain and a Davidson potential which reduces sp(3,R) D [R^]so(3) Z) so(3). Thus, we expect 
results that lie between the two limits. In fact, the results of our calculations prove to be such that both [M^]so(3) 
and su(3) are extremely good quasi-dynamical symmetries. This is a remarkable result, because although su(3) and 
the rigid-rotor algebra have similar algebraic structures, they have very different physical interpretations. 

A Lie algebra g is said to be a quasi-dynamical symmetry for a Hamiltonian H if the matrix elements of q between 
a set of eigenstates of H are equal to those of an irrep of g, even though the eigenstates in question do not belong to 
a subspace of the Hilbert space for H that is irreducible under the action of g. At first sight this would be appear to 
be an unlikely physical situation. Indeed, if matrix elements of an algebra g were found in some situation to be equal 
to those of an irrep, one might be tempted to infer that g is a full dynamical symmetry for the Hamiltonian. The soft 
rotor model is the prototype of a situation where this is not the case, as we now show. 

Let {^KLMi(3,j)} denote a basis of wave functions for a rigid rotor with intrinsic quadrupole moments 

Qo^kPcosj, Qi = Q_i=0, Q2 ^ Q-2 ^ —i=kl3ain-f , (3) 
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where fc is a suitable constant, cf. Eq. (g^. Then soft-rotor wave functions can be expressed 

'i'aKLM = / /a (/3, 7) '^KLM (/?, 7) dw , (4) 



where dv is a suitable volume element for /? and 7. One finds that the matrix elements between all states of the same 
a are equal to those of a rigid rotor of deformation 

(Oo)=y"fc/3cos7/„(/3,7)d«, {Q2)^ J ^kl3 sin J f4f3,j)dv. (5) 

Thus, the states of a given a span an embedded representation of the rigid-rotor algebra in the terminology of ref. |17| ] 
and we say that the rigid-rotor algebra is a quasi-dynamical symmetry for the soft rotor. 

An examination of the physics of the situation quickly reveals that the rigid-rotor algebra is a quasi-dynamical 
symmetry whenever the rotations are adiabatic relative to the complementary vibrational motions. Thus, one obtains 
the rigid-rotor algebra as a quasi-dynamical symmetry for a rotor-vibrator Hamiltonian whenever the rotational- 
vibrational (i.e., Coriolis and centrifugal) coupling interactions are negligible or omitted. They can be taken into 



account subsequently, as quasi-dynamical symmetry-breaking perturbations, as is standard in the pheomenological 
rotor model. 

Now it is known [|8[ that the su(3) algebra contracts to the rigid-rotor algebra [IR^]0o(3) and that this contraction 
is realized for large-dimensional irreps of su(3). Thus, we say that su(3) is a quasi-dynamical symmetry for a sequence 
of states if their expansions on an su(3) basis are of the form 

*c./fLAf = ^aAM*KV> (6) 

A/i 

where {^j^j^m} ^^ ^ basis of states of a large-dimensional su(3) irrep of highest weight (A, /i). Note that shell-model 
states can always be expanded on an su(3) basis. However, we only obtain su(3) as a quasi-dynamical symmetry if 
the coefficients {Cax^i} are independent of KLM for a useful range of values of the latter. 

The fact that su(3) is a good quasi-dynamical symmetry for major-shell mixing Hamiltonians is extremely important 
for achieving the eventual goal of including sp(3,R) symmetry-breaking interactions in the symplectic model. For 
although, interactions like the spin-orbit and pairing interactions break sp(3,M), they have been shown to preserve 
quasi-dynamical su(3) symmetry ||5|,|7[|. This is expected to occur, to a good approximation, whenever the rotational 
motions one is describing are known from experiment to be adiabatic. 

The concept of quasi-dynamical symmetry can be regarded simply as a group-theoretical expression of the standard 
methods for handling adiabatic decoupling of collective motions along the lines of the Born-Oppenheimer approxi- 
mation |20(] and the standard nuclear rotor model. A formulation in the precise language of group theory has the 
advantage that it opens the concept up to more general application. An overview of the concept was given in a recent 
conference report js). 

In Sect. O, the Hamiltonian of the nuclear symplectic model is described and the phy sical motivation behind the 



model in connection with Bohr-Mottelson collective model Il2| is explained. Sect. HI defines the sp(3,]R) and its 



subalgebras relevant to this analysis. Sect. [V describes the construction of basis states and matrix elements within 
the three dimensional harmonic oscillator space. The results are given in Sect. M and some conclusions are drawn in 



Sect. V] 



II. THE SYMPLECTIC MODEL HAMILTONIAN 

The Hamiltonian of the nuclear symplectic model, in its simplest form, consists of two parts: 

H = Ho + ViQ) ; (7) 

a three dimensional many-particle harmonic oscillator (shell-model) Hamiltonian 






(8) 



and a collective potential V{Q), where Q is the quadrupole tensor for the nucleus. 

The collective potential V{Q) is a rotationally-invariant function of the quadrupole moments. In a Cartesian basis, 
the quadrupole moments for a nucleus are given by 

Q^J = Qu- - ^S,, ELi Qkk , (9) 

where Q is a monopole-quadrupole tensor with components 

A 



A A 



n— 1 m,n— 1 



and {Xi] are the components of the center-of-mass vector X = ^ ^^ x„. Removal of the center-of-mass contribution 
to the quadrupole moments in this way, ensures that the spurious A^-phonon center-of-mass states of the Hamiltonian 
remain, unmixed with other states, at an excitation energy of Nfuo. 



A general, rotationally-invariant potential can be expressed 

V{Q) = V{Q-Q,Q-QxQ), (11) 

as a function of the quadratic and cubic scalars 



2 



Q-Q= }] (-I)^Q-mQm « TrQ 



where 



2 

Q-QxQ= Yl (-1)^(21^; 2p|2/^)Q_^Q,Qp ex TrQ3, (12) 



Qt^=Yl ^n^2M(i-n) 



with r„ = |x„ — X| and f„ = (x„ — X)/r„. In the study reported in this article, the collective potential is taken to 
be the Davidson potential 

yiQ)-x{Q-Q + ^). (13) 

This potential is shown as a function of Qq — y/Q ■ Q in Fig. ^. 

The value of the parameter e determines the value of Qo ^^ which the potential has its minimum value. Thus, 
e is chosen such the potential has a minimum at the observed deformation of the nucleus under investigation. The 
strength x of the potential is then set such that the lowest energy {L — 0) wave function that emerges is such that 
the expectation value {Q ■ Q) in this state is equal to the value of Qq at which the potential is a minimum. We refer 
to this as the self-consistent value of x- 

III. THE SYMPLECTIC ALGEBRA AND ITS SUBALGEBRAS 

The quadrupole moments {Q^i}, the many-nucleon kinetic energy J2nPn/'^''^j ^^'^ the harmonic oscillator poten- 
tional ^mui^ ^„ r^ are all elements of an sp(3, R) Lie algebra. Thus, sp(3, M) is a spectrum generating algebra for H; 
i.e., it is the Lie algebra of a dynamical group for H. It folows that the eigenstates of H belong to a single irreducible 
representation of Sp(3, E). 

The complex extension spc(3,IR) of sp(3,R) is spanned by the operators (in a Cartesian coordinate system) 

A A 

n— 1 m.n— 1 



^ij — / ^ ^7ii^7ij A / ^ ^7ni^7ij 1 V 7 

n— 1 m.n— 1 

^ A ^ A 

71—1 771, n — 1 

where b]^- and bni are the diniensionless harmonic oscillator raising and lowering (boson) operators 

, _J_f J_ \ 

^7ii , /—I ^ni + Pni I ; 

5ov2 V mu} ) 

in units of the oscillator length fep = {h/mujY''^. The latter operators satisfy the commutation relations 



[&„,I] = [foL,II]-0, (15) 

of a Heisenberg-Weyl algebra. 

The symplectic algebra sp(3,R) contains many subalgebras, including the u(3) D su(3) chain of Elliott's SU(3) 
model 11^ and the rot(3) {— [R^]5o{3)) rigid-rotor algebra of Ui's model jl^]. The u(3) subalgebra is spanned by the 
{Cij} operators. One sees, for example, that u(3) contains the harmonic oscillator Hamiltonian Hq = huj'^^^^ C'u as 
an element. The rot(3) algebra is spanned by three angular momentum operators 

Li = ~i{C.jk - Ckj) , (i, j, k cyclic) , 

and five components {Q^; IJ. — 0, ±1, ±2} of the L — 2 quadrupole tensor Q. 
The Cartesian components of the monopolc-quadrupole tensor are expressed 

Qy = bl [Ay + i (Q, + C,,,) + B,j] . (16) 

The spherical components of the quadrupole operators are 

Q^ = blV3{A2^ + C2t,+B2^). (17) 

One sees that the part of Q^ that commutes with Hq is the su(3) quadrupole operator 

Q^ = blV3C2^. . 

Thus, Bu(3) can be viewed as the projection of the rot(3) algebra onto the space of operators that leave spherical 
harmonic oscillator shells invariant. 

An irrep of u(3) is characterized by a highest weight state \(p) and a corresponding highest weight uj — {uji,uj2,0J3) 
defined such that 

C^j\(|)) ^ ior i < j , C^^\(f>) ^ UJ^\(j)) . 

Such an irrep remains irreducible on restriction to its u(l) + su(3) C u(3) subalgebra and has u(l) + su(3) highest 
weight iV(A^), where 

N ^ UJi + LU2 + ^3 , X ~ UJi ~ UJ2 , /i = aJ2 - tJ3 . 

An irrep of the rigid-rotor algebra, is characterized by an intrinsic state which is an eigenstate of the quadrupole- 
moment operators with eigenvalues {Qfj,} that are related to the shape variables of the Bohr-Mottelson collective 
model |l|] by 

^0 = v/^^-Ro/3cos7, 



Q±2 = ^^ARlP^^^l, (18) 

Q±i = , 

where A is the mass number, Rq is a nuclear radius, and (3 and 7 are deformation and asymmetry parameters, 
respectively. 

It is useful to note that the elements of sp(;;(3,M) are all compnents of U(3) tensors, i.e., they transform according 
to irreducible representations of u(3); the operators {Ay} are components of an irreducible tensor A of highest weight 
{200}, the {Bij} are components of an irrreducible tensor B of highest weight {00 — 2}, the operators {Cij , Cji.Ca — 
Cjj'ji < j} are components of an irreducible su(3) tensor C'^^^^ of su(3) highest weight (11) (u(3) highest weight 

{10 — 1}), and Hq — fku J2i=i ^a i^ ^ u(3) scalar; a tensor of highest weight {000}. 

Now observe that the first term, Hq, of the Hamiltonian of Eq. (R) is SU(3) (and U(3)) invariant whereas the second 
term is invariant under the dynamical group R0T(3) {— [R^]S0(3)) of a rigid rotor. Thus, the two components, Hq 
and V{Q), of the Hamiltonian respectively reduce the two subalgebra chains: 



sp(3.R) 

/ \ 

su(3) rot(3) 

\ / 

so(3) 

This implies that the eigenstates oi H = Hq + V{Q) are intermediate between those of the SU(3) and rigid-rotor 
models. We show in the following that, in fact, both SU(3) and R0T(3) are remarkable good quasi-dynamical 
symmetries for this Hamiltonian. 

IV. BASIS STATES AND MATRIX ELEMENTS 

A unitary irrep of sp(3,R), within the shell- model space of a mass- A nucleus, is characterized by a lowest- weight 
state |crlw), with weight a = (cri, (72, (73), defined by the equations 

Bij\a\w) = , Cy |crlw) = for i < j, Cii\a\w) = cr,;|crlw) . (19) 

When (7 is a triple of positive integers or of positive half-odd integers, the corresponding irrep is either a discrete 
series representation or, for mass number ^ < 6, a limit of a discrete series irrep. 

Let {|cra)} denote an orthonormal basis for the subspace of states of an sp(3,M) irrep, of lowest weight a, that 
satisfy the equation 

By|cra)=0. (20) 

We refer to these states as vacuum states for the corresponding sp(3,]R) irrep. They are a basis for a u(3) irrep of 
highest weight a. Moreover, it is known that a basis for an sp(3, M) irrep can be constructed by acting on the vacuum 
states with polynomials in the {^y} raising operators. Let 

41m(^) = [^x^X---X^]^"1m (21) 

denote a product of N symplectic raising operators coupled to the {KLM) component of a U(3) tensor operator 
Z(")(v4) of highest weight n — [rii, n2,'n^), where rii, 712, and 713 run over the even-integer values for which 

ni>n2>n2>0, ni -f n2 -f n3 = 27V . (22) 

Acting on the vacuum states with these tensor operators gives basis states 

\anpojKLM) = [Z^^^A) x \a)]''^^^, (23) 

which reduce the subgroup chain 

Sp(3,M) D U(3) D S0(3) D S0(2) 

a np u K L M ' ^ ' 

These states are eigenstates of the harmonic oscillator hamiltonian 

Ho\{cTn)pujKLM) = {Nq + ni+n2 + n3)huj\{an)pujKLM) , (25) 



where Nq — cti + a2 + (T3. However, although the basis states of Eq. ( |23D span the carrier space of an sp(3,M) irrep, 
they are not an orthonormal basis due to a multiplicity of U(3) subirreps labelled by the indices n and p. Thus, it is 
necessary to take suitable linear combinations 

larcuKLM) = Y,iK)nlr\<^npujKLM) 

np 

to form an orthonormal basis. The K matrix coefhcients are conveniently determined by vector coherent state (VCS) 
methods |l^ , p5[ . The calculation of matrix elements of elements of the sp(3,M) Lie algebra in the corresponding 
orthonormal basis is also straightforward using VCS theory. 

Matrix elements of the components {Xkim} of an SU(3) tensor X^^"^' are obtained from their su(3)-reduced matrix 
elements using the generalized Wigner-Eckart theorem 



{av' t' uj' K' L' M'\Xkira\(JTUjKLM) 



7 



{{Xfi)KL; {pq)kl\\j{X'fi')K'L') {LM; lm\L'M') 



x(crTV||X(P9)| 



11 ■> 



(26) 



where X = uji — W2, fi — 102 — W3, and 7 indexes the niultiphcity of SU(3) irreps of highest weight {X' 11') in the tensor 
product (A^) X (pq); {LM;lm\L'M') is an S0(3) Clebsch-Gordan coefficient and {{X^)KL; {pq)kl\\^{X' ^')K' L') is an 
S0(3)-reduced Clebsch-Gordan coefficient for SU(3). 

SU(3)-reduced matrix elements of elements of the sp(3,M) Lie algebra are given, for example, in ref. jlj]. To 
calculate matrix elements oi Q ■ Q, following Rosensteel fell, we start by normal ordering the expansion 



= 3(^2 + C2 


+ B2) ■ iA2 


+ C2 


+ B2 


) 




= 6Cf (^^ - 


3L^ + lOHo + IOVQBo 


+ 6A2 


■B2 + 


ion relations 


•> 












B2 


■A2- 


-A2 


B2 = 


10 TT 




B2 


■C2- 


-C2 


S2- 


fV6B 



{6C2 ■ B2 + W2 ■ B2 + h.c.} , 



• 



(27) 



(28) 



With such a normal-ordered expansion, i.e., with B2 operators on the right and A2 operators on the left, we do 
not have to include intermediate states external to the truncated space; consequently, the intermediate sums are 
minimized and the calculations are less time consuming. Optimization of the computations in this way is important 
because the number of basis states grows exponentially as the number of major oscillator shells in the calculation 
increases. For example, suppose we want to calculate matrix elements of Q • Q between states of a truncated space 
comprising the oscillator shells Nohuj, {Nq + 2)tiuj, {Nq + 4:)huj, . . ., {Nq + n)Huj. Without normal ordering, we would 
have to include intermediate states in the calculation from the (A^o + n + 2)Hlu shell. The number of extra states 
involved in the calculation would then be proportional to n(n -t- 1). 

The first three terms in Eq. (p7|) are diagonal in the chosen basis with eigenvalues given by 



5U(3)\ 



E^i^—kN){u;, 



iiV- 

3 



2i) = f (A^ +fi^ + Xfi + 3X + 3/i) , 



{L^)=L{L + 1), {Ho)=N. 

The quadratic terms are expressible as su(3)-coupled tensors using the identity 



p(Am) 



p(A/i) 




(29) 



(30) 



this gives 



A2-B2 



IVeiA X i3)("0) + iV30(A X B)('') , 



,(02) 



C2 • B2 = V5(C X B 
B2-B2 = f V5(B X B)^' 



(31) 



(04) 



liBxB)^^'l 



The su(3)-scalar component of A2 ■ B2 is related to the quadratic Casimir invariants of sp(3,M) and su(3) 



(AxB) 



(00) 



^V6{C: 



.SU(3) 



1 tt2 



4iJn - C: 



.Sp(3J 



'). 



where 



with 



^^sp(3,K)^ ^ J2^ ^^(^^ _ 2t) = |(A2 + AoMo +f^l + 3Ao + 3mo) + ^N^ - ANo , 



Xo — Cri — (72 , Mo = 0-2 

Combining the above results, Eq. (E7h becomes 



(T3 , Nq = (Ji + a2 + 0-3 . 



(32) 



(33) 



g.Q = -|cf(^-* 



17^SU(3) 07-2 



6^0 



30(^ X 5)1,^^^ 



-{IOV6B0 + 6%/5(C X B)J,"^^ 
-2V5(B X i?)(°*) + 5{B X bC^ 



h.c.}. 



(34) 



V. RESULTS 

Calculations have been done for a typical heavy nucleus with oscillator energy huj — 7.49 MeV; a value appropriate 
for the ^|g Er nucleus. The deformation parameter e — 1.01 x 10^^ gives the minimum of the collective potential at 
/3in — 0.35; a value close to that inferred from the experimental B{E2:2f -^ 0^) transition rate. We choose an sp(3, M) 
irrep with lowest weight (327^,249^,249^) (or equivalently, 826^(78,0) in a U(l) x SU(3) notation §). This irrep 
is deduced from an empirical formula for the intrinsic mass quadrupole moment of a deformed oscillator |22|Jl5|| . The 



strength x of the potential V{Q) in Eq. (13) is varied in this study. The ratio between the potential and oscillator 
strengths determines the structure of the wave functions. Due to computational limitations, the space is restricted to 
states belonging to shells below 12hLu; apart from this restriction there is no further truncation. 

A. Ground state band 

Energy spectra for different x values (in units of MeV) are given in Fig. y. The results show ground state bands 
with rotational spectra and moments of inertia that decrease as x increases. Moments of inertia, defined for each 
angular momentum state, by the expression 



'^-^slir^^^^^ + i) (^^) 



are shown in Table |[ 



TABLE L Moments of inertia for states of angular momentum L; the values in parenthesis are percentages relative to the 
value at L = 2. 



L\X 


1.67x10"^ MeV 


3.33x10-*^ MeV 


1.67x10-* MeV 


8.33x10"'' MeV 


2 


207.62 MeV-^ (100.00%) 


148.27 MeV"' (100.00%) 


70.10 MeV-^ (100.00%) 


15.72 MeV-^ (100.00%) 


4 


207.44 (99.91%) 


148.15 (99.92%) 


70.04 (99.91%) 


15.70 (99.89%) 


6 


207.16 (99.78%) 


147.97 (99.80%) 


69.94 (99.76%) 


15.68 (99.72%) 


8 


206.78 (99.59%) 


147.72 (99.63%) 


69.79 (99.55%) 


15.64 (99.49%) 


10 


206.29 (99.36%) 


147.41 (99.42%) 


69.61 (99.30%) 


15.59 (99.20%) 


12 


205.71 (99.08%) 


147.02 (99.16%) 


69.40 (98.99%) 


15.54 (98.85%) 


14 


205.02 (98.75%) 


146.57 (98.86%) 


69.14 (98.62%) 


15.47 (98.43%) 


16 


204.22 (98.36%) 


146.06 (98.51%) 


68.84 (98.20%) 


15.40 (97.95%) 


18 


203.33 (97.93%) 


145.48 (98.12%) 


68.50 (97.72%) 


15.31 (97.41%) 



The most remarkable result is that, for each of the interaction strengths shown, the spectra are almost identical 
to those of a rigid rotor with excitation energies very accurately proportional to L{L + 1). B{E2) transition rates 
between adjacent states, shown for x = 3.33 x 10"* MeV in Table |l[ are also in remarkably close agreement with 
those of a rigid rotor. These results are significant because, they give rotor model results with fully microscopic 
166-particle wave functions. Thus, they provide us with the means to explore the rotational dynamics of a nucleus 
at the microscopic, many-nucleon, level. As a comparison, the results from Elliott's SU(3) model are also given with 
the effective charge e = 2.04. In the symplectic-Davidson model, e = 1. 

A question of considerable interest is the nature of nuclear rotational energies. In the Bohr-Mottelson rotor model 
| jl2[ rotational energies are interpreted as arising from the kinetic energy whereas in Elliott's SU(3) model [Q they 
come from the potential energy. Since the symplectic model contains both a rigid-rotor and Elliott's SU(3) models as 
limiting submodels, it is of considerable interest to see how it interpolates between the two limits. In the symplectic 
model, the kinetic energy operator is written 



Ho~JUAo + Bo) . (36) 



The contribution of the kinetic to the total calculated excitation energy is shown as a percentage for each state of the 



ground state band in Table Q). It can be seen that, for the smaller values of X: the kinetic energy gives a negative 
contribution to excitation energies. Its contribution is positive for larger values of x, but remains much less than that 
of the potential energy. Note, however, that because of the truncation to shells of 12fiuj and below, the results shown 
for X = 8.33 X 10"^ MeV are not fully converged and are unreliable. Even so, since we obtain rotational bands in 
close agreement with the rotor model for a wide range of potential strengths, it is hard to avoid the inference that 



TABLE II. The reduced quadrupole transition strength B{E2 : L -^ L—2) in Weisskopf unit (W.u.) for the experimental 
data, SU(3) model, this model (with coupling constant x = 3.33 x 10^* MeV), and a rigid rotor. Ratios of transition strengths 
to those of the the 2^0 transition are given in parenthesis. 



Li -^ Lj 


experiment 


SU(3) 


Davidson 


rigid rotor 


2^0 


214 ± 10 (1.00) 


207 (1.00) 


207 (1.00) 


207 (1.00) 


4-^2 


311 ±20 (1.45) 


295 (1.43) 


295 (1.43) 


296 (1.43) 


6->4 


347 ±45 (1.62) 


324 (1.57) 


325 (1.57) 


326 (1.57) 


8^6 


365 ± 50 (1.70) 


338 (1.63) 


339 (1.64) 


340 (1.64) 


10^8 


371 ±46 (1.73) 


345 (1.68) 


348 (1.68) 


350 (1.69) 


12^ 10 


376 ± 40 (1.76) 


349 (1.69) 


353 (1.71) 


356 (1.72) 


14^ 12 




350 (1.69) 


356 (1.72) 


361 (1.74) 


16 -> 14 




350 (1.69) 


357 (1.73) 


365 (1.76) 


18^16 




349 (1.69) 


358 (1.73) 


367 (1.77) 



TABLE III. The kinetic energy as a percentage of the total excitation energy for different values of the coupling constant x. 
The values of /3o show the average deformation of the (L = 0) ground state according to Eq. (|lq). 



L\X 



1.67x10"^ MeV 



3.33x10^" MeV 



1.67xlO~''MeV' 



8.33x10"'' MeV 



2 
4 
6 
8 
10 
12 
14 
16 
18 



-31.95% 
-31.93 
-31.90 
-31.85 
-31.79 
-31.72 
-31.64 
-31.54 
-31.43 



-21.57% 
-21.56 
-21.53 
-21.49 
-21.44 
-21.37 
-21.30 
-21.22 
-21.12 



26.43% 
26.44 
26.46 
26.48 
26.52 
26.55 
26.60 
26.65 
26.71 



9.46% 
9.46 
9.45 
9.44 
9.43 
9.41 
9.39 
9.37 
9.35 



/3o 



0.247 



0.275 



0.345 



0.351 



nuclear rotational energies are most likely not 100% kinetic in origin. If this is correct, it has considerable conceptual 
implications for the interpretation of nuclear rotational dynamics. We may continue, for convenience, to describe the 
cofactor of the L{L + 1) rotational energy as the inverse of a "moment of inertia", but one should recognize that, if 
the excitation energy is not kinetic, the concept is misleading. 

Another result of this study is that the moments of inertia as well as the kinetic energy portions of the excitation 
energy change less than 5% over the range of angular-momentum values considered. This result can be interpreted 
as signifying that the states of a band have a common intrinsic structure that changes little with increasing angular 
momentum. This interpretation becomes much more compelling when one observes the behaviour of the coefficients 
of the wave functions in the expansion 



(Am) 



(37) 



of the states of the ground state band in a U(3) Z) SU(3) D S0(3) basis. The coefficients are plotted for different x 
values in Fig. |^. The figure shows clearly, that for all values of x considered, the coefficients are essentially independent 
of L and negligible for K ^ 0. Moreover, for the larger %, this is in spite of a huge mixing of SU(3) irreps from many 
major shells. It follows that, while SU(3) is far from being a good dynamical symmetry, it remains an extrordinarily 
good quasi-dynamical symmetry according to the definition given in the Introduction and Refs. [[7|Jl7|. 

Observe also that the distribution of U(3) irreps is dominated by the so-called stretched irreps; the stretched irreps 
are those of the sequence 



A^o(Ao,0), iVo ± 2(Ao ± 2, 0) , TVq ± 4(Ao + 4, 0) ^ 



iVo + 2n(Ao + 2n, 0) , 



The value of the coupling constant for which the L = ground state of the model Hamiltonian has a deformation 
/So = 0.350 equal to /?in, the value for which the Davidson potential is a minimum, was found by repeated calculation 
to be given by x = 3.33 x 10"'* MeV. We call this the self-consistent coupling constant. The spectrum of the ground 
band for this x is displayed in Fig. H in comparison with the observed spectrum of ^^^Er and that of the rigid rotor. 
One sees that the results track the rigid rotor more closely than they do experiment. This may be an artifact of 



the Davidson potential which continues to rise unrealistically at large deformation and excessively inhibits centrifugal 
stretching. We also computed the reduced electromagnetic transition strengths 

The ratios of transition strengths, shown in parenthesis in Table Q and also depicted in Fig. 0, are in extremely 
close agreement both with experiment and the rigid-rotor model. Their magnitudes shown in Fig. || are also in good 
agreement with experiment. 

B. Giant resonance bands 

In addition to the ground-state rotational band, the symplectic model gives many excited bands. However, for a 
single iVo(Ao,0) irrep, as considered here, without spin-orbit, pairing and other irrep-mixing interactions, the simple 
symplectic model has no low-lying excited bands any more than the SU(3) model has excited bands for a (Aq, 0) irrep. 
The lowest-energy excited bands of the simple symplectic model are associated with the giant monopolc (breathing 
mode) and giant quadrupole resonance degrees of freedom. 

For small values of X: these occur in the model, with a Davidson potential, at around 2huj. They are shown for 
several values of x in Fig- H; results for the self-consistent value are given in Fig. pi 

Two results are worth noting. The first is that the energies of the giant-resonance bands rise with increasing values 
of X and become unrealistically high at the value considered appropriate for the ground state band. This we believe to 
be a reflection of the fact that, although the Davidson potential has many useful features, it rises too steeply away from 
the equilibrium deformation. The second notable result is that the symplectic model gives three monopole-quadrupole 
giant resonance bands; two of these, the K ^ and K ~ 2 bands, occur also (albeit at much lower energy) in the 
Bohr-Mottelson model where they are associated with beta- and gamma-vibrations, respectively. However, whereas 
there is no iiT = 1 one-phonon band in the Bohr-Mottelson model, such a band is non-spurious in the symplectic 
model which includes intrinsic vorticity degrees of freedom. 

The monopole and quadrupole giant resonance states in real nuclei are believed to lie at energies close to 2tiuj. 
However, their strength is invariably fragmented and much of it lies in the continuum. Thus, we make no attempt to 
compare our results for the GR states with experiment. However, it is of interest to examine the structure of the wave 
functions. Fig. ^ shows the decomposition of the K = (giant-beta band) states in terms of their SU(3) components, 
for the self-consistent value of x, in comparison with the corresponding decomposition of the ground-band states. It 
can be seen that the amplitude coefficients (including their signs) are independent of angular momentum to a high 
degree of accuracy. It can also be seen that the coefficients for the two sets of states have the same signs except for 
those of the stretched states, which are of opposite sign. Note that the relative signs of the amplitude coefficients for 
different SU(3) subirreps are determined by SU(3) Clebsch-Gordan coeffients 124| and have no particular meaning, 
in general, except for the stretched states which all have the same sign. However, a change in the relative signs 
of coefficients between the ground and excited states is meaningful. Thus, it is fortuitous that the coefficients for 
the stretched states all have positive sign for the states of the ground state band because it highlights the fact that 
these coefficients change sign for the giant beta band. On reflection this is what one would expect for a giant-beta 
vibrational excitation (see Fig. |8|). 

VI. CONCLUSION 

The symplectic model is currently the only model that is capable of giving rotational states for heavy nuclei as 
eigenstates of a rotationally-invariant Hamiltonian in a realistic shell-model space. Thus, the ability of the model to 
give the energy levels and E2 transitions strengths between states for ground-state rotational bands, without the use 
of an effective charge, provides a powerful framework for understanding the dynamics of nuclear rotations in terms of 
interacting neutrons and protons. 

Early applications of the model by Park et al. [^ were remarkably successful and raised many interesting questions. 
However, because of the severe truncations of the space that were necessary at the time, the reliability of the results 
could be questioned. In particular, one was concerned that the ability of the model to give correct moments of inertia 
might be lost on increasing the size of the model space. The results of the present calculations show that this is not 
the case. 

A particularly interesting challenge was to learn how a model, without pair correlations, could give correct moments 
of inertia when it is known that the cranking model is only succesful when pairing correlations are included. The early 

10 



calculations of Park et al. indicated that the dominant contribution to rotational energies came from the potential 
energy part of the Hamiltonian, thus calling into question the very concept of the moment of inertia as an inverse 
coefficient of the L^ term in the kinetic energy. The results of the present calculation indicate that the inclusion of 
only stretched states, as in the calculation of Park et al., tends to exagerate this effect. Nevertheless, it confirms that 
the dominant component of the rotational energies comes from the potential energy; for the self-consistent value of x 
only about 20% of the rotational energy comes from the kinetic energy in the present calculation. 

In addition to investigating rotational states up to much higher angular momentum {L — 18), the present calcu- 
lations have focussed on understanding the structure of rotational states in terms of their SU(3) content. We have 
shown that although there is huge mixing of SU(3) irreps from many major harmonic oscillator shells, the mixing 
is highly coherent and establishes SU(3) as a remarkably good quasi-dynamical symmetry for the model. Other cal- 
culations li~ 



reviewed in [pi, show that such quasi-dynamical symmetry is also conserved when SU(3) irreps are 
further mixed by spin-orbit and pairing forces. Thus, the results show that, as far as the calculation of E2 transition 
rates are concerned, the use of a single SU(3) irrep with an effective charge will give accurate results. However, for 
other observables, not related to elements in the SU(3) algebra, there is no reason whatever to expect effective charge 
methods to take account of the large mixing of SU(3) irreps observed. 

The close agreement between the results of the symplectic model calculation and those of the rigid-rotor model, 
shows that, by definition, the rigid-rotor algebra is also an excellent quasi-dynamical symmetry for the symplectic 
model with a Davidson potential. But again, for observables not related to the rigid-rotor algebra, such as electron 
scattering current operators, one cannot predict what the results will be by purely algebraic methods. 

The fact that two competing dynamical subgroup chains, although very different in their physical content, can both 
be good quasi-dynamical symmetries is remarkable but understandable. It can happen because, for large-dimensional 
representations and states of relatively low angular momentum, an SU(3) irrep contracts to an irrep of the rigid-rotor 
algebra |Lq |. Thus, it transpires that the lower angular-momentum states of a large-dimensional SU(3) irrep belong 
to an embedded representation |1^ of the rigid-rotor algebra and vice-versa. 

Having demonstrated that symplectic model calculations with phenomenological (albeit microscopically expressable) 
potentials have the ability to describe nuclear rotational states, our next goal would be ideally to perform calculations 
within the same (large) shell model space but with realistic two-nucleon interactions. Such calculations can and 



have been contemplated for light nuclei |19|. But they are computer intensive and impractical for heavy rotational 



nuclei. We would even like to go further and include spin-orbit and short-range (e.g., pairing) interactions which 



mix different sp(3,M) irreps. We would like to carry out calculations for superdeformed bands 1 26 2^ as well as for 
normally deformed low-lying rotational bands of heavy nuclei. Superdeformed bands are naturally associated with 
excited representations of the symplectic model which fall into the low-energy domain (as do Nilsson model states) 
as a consequence of shell effects in a deformed shell model. The challenge is to explain why they do not mix more 
readily with the large density of less deformed low-energy states. We believe that symplectic symmetry and su(3) 
quasi-dynamical symmetry has the potential to answer such questions. 

It is unlikely that such calculations will ever be done in a shell-model space sufficiently large to ensure convergence 
of the results. However, the results of the present study and previous investigations of mixing SU(3) irreps with 
spin-orbit |17[ and pairing interactions M imply that, when experiment finds states of a nucleus that are fitted well by 
the rotor model and to have a large deformation, then we have reason to believe that SU(3) is a good quasi-dynamical 
symmetry for these states. Armed with this information, one can hope to design realistic mixed SU(3) calculations 
within large spaces. In particular, one can calculate just one representative angular momentum state for each band, 
with the understanding that the coefficients should be the same (to a good approximation) for all other states of 
the band. Alternatively, one could carry out calculations in an intrinsic space which includes just the highest weight 
states for the contributing SU(3) irreps. 
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FIG. 1. The Davidson potential; its two components are shown separately. 
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FIG. 2. The energy spectra for various strengths of the potential. The dashed hnes connect energy levels of the same angular 
momenta. 
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FIG. 4. Energy spectra for experimental data, the SU(3), the symplectic-Davidson, and the rigid rotor models. The B(E2) 
values are given inside the boxes. 



16 



2.0 



^1.5 
o 

A 

I 

c\i 

4i 
CD 

^1.0 

_i 

A 

I 



(M 

LU 



CD 



0.5 



0.0 




-♦--experiment 
■A" SU(3) 
-0^ Davidson 
-▼- rigid rotor 



12 



16 



20 



FIG. 5. Relative B{E2) values for experimental data, the SU(3), the symplectic-Davidson, and the rigid rotor models. 



17 



keV 



40000 



30000 



20000 



10000 



5 hco 



K = (GMRj 



K = 2 (GQR) 



K = 1 



ground 
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FIG. 8. As for Fig. |7| except that only L — coefficients are plotted and the phase factor for the ground-state coefficients are 
set to be positive. The stretched bases are indicated by asterisks. The insets show changes of phase for the first few coefficients. 
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